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Abstract. The ep -^ e'pp^ reaction has been measured, using the 5.754 GeV electron beam of Jefferson Lab 
and the CLAS detector. This represents the largest ever set of data for this reaction in the valence region. 
Integrated and differential cross sections are presented. The W, Q^ and t dependences of the cross section 
are compared to theoretical calculations based on i-channel meson- exchange Regge theory on the one hand 
and on quark handbag diagrams related to Generalized Parton Distributions (GPDs) on the other hand. 
The Regge approach can describe at the ~ 30% level most of the features of the present data while the 
two GPD calculations that are presented in this article which succesfuUy reproduce the high energy data 
strongly underestimate the present data. The question is then raised whether this discrepancy originates 
from an incomplete or inexact way of modelling the GPDs or the associated hard scattering amplitude 
or whether the GPD formalism is simply inapplicable in this region due to higher-twists contributions, 
incalculable at present. 

PACS. 13.60.Le Production of mesons by photons and leptons - 12.40.Nn Regge theory - 12.38.Bx 
Perturbative calculations 
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1 Introduction 

The exclusive electroproduction of photons and mesons 
on the nucleon is an important tool to better understand 
nucleon structure and, more generally, the transition be- 
tween the low energy hadronic and high energy partonic 
domains of the Quantum Chromodynamics (QCD) theory. 
Among all such exclusive processes, the ep -^ e'pp^ 
reaction bears some particular advantages. It is a pro- 
cess for which numerical calculations and predictions are 
available both in terms of hadronic degrees of freedom, 
via Reggeized meson exchanges, and in terms of partonic 
degrees of freedom, via Generalized Parton Distributions 



(GPDs). We refer the reader to refs. [I||2|[3] and refs. [H 
[51[M7[|^[|^[TU] for the original articles and general reviews 
of Regge theory and GPDs respectively. Defining Q^ as 
the absolute value of the squared mass of the virtual pho- 
ton that is exchanged between the electron and the target 
nucleon, partonic descriptions are expected to be valid at 
large Q^, while hadronic descriptions dominate in photo- 
and low-Q^ electroproduction. Fig. [1] illustrates these two 
approaches at the electron beam energies available at Jef- 
ferson Laboratory (JLab). Concerning the Reggeized me- 
son exchange approach, the total and differential cross 
sections associated with the exchanges of the dominant 



Regge a and /2 trajectories have been calculated by Laget 
et al. [TT|ll2j. Concerning the GPDs approach, the so- 
called "handbag" diagram, with recent modelings of the 
unpolarized GPDs, has been calculated by two groups: 
Goloskokov-Kroll [13] and Vanderhaeghen et al. P^irMTOl 
[TT] . Let us note here that in the GPD approach the leading 
twist handbag calculation is valid only for the longitudi- 
nal part of the cross section and that, experimentally, it 
is important to separate the longitudinal and transverse 
parts of the cross sections when measuring this process. 



7l,t 




Fig. 1. The mechanisms for p" electroproduction at JLab ener- 
gies for low Q^ (left diagram) through the exchange of mesons 
and for high Q^ (right diagram) through the quark exchange 
"handbag" mechanism (valid for longitudinal photons) where 
H and E are the unpolarized GPDs. 



This article presents results for the exclusive electro- 
production of the p^ vector meson on the proton mea- 
sured with the 5.754 GeV electron beam of the CEBAF 
accelerator and the CEBAF Large Acceptance Spectrom- 
eter (CLAS) at JLab. The aim of this analysis is to com- 
pare the integrated and differential cross sections of the 
ep —>■ e'pp^ reaction that have been extracted over the 
intermediate Q^ region accessible at CLAS, with the two 
Regge and GPD theoretical approaches, and thus deter- 
mine their domain of validity and constrain their various 
inputs. 

There are a few existing electroproduction data in a 
similar kinematical regime: early data with the 7.2 GeV 
beam at DESY ^ and with the 11.5 GeV beam at Cor- 
nell [in], and more recently with the 27 GeV beam at 
HERMES piT and the 4.2 GeV beam of JLab [H]. The 
present work explores new phase-space regions and, in re- 
gions of overlap, has much finer binning and precision. 

In section [2] we present the experimental procedure we 
have adopted to extract our integrated and differential 
cross sections. In section [3l after briefly describing the 
Regge and GPDs models, we compare these calculations 
to our data. Finally, we draw our conclusions in section |4j 



2 Experimental procedure 

The CLAS detector ^ is built around six superconduc- 
tion coils that generate a toroidal magnetic field primarily 
in the azimuthal direction. Each sector is equipped with 



three regions of multi-wire drift chambers (DC) and time- 
of-flight scintillator counters (SC) that cover the angular 
range from 8° to 143°. In the forward region (8°< 9 < 
45°), each sector is furthermore equipped with a gas- filled 
threshold Cerenkov counter (CC) and a lead-scintillator 
sandwich type electromagnetic calorimeter (EC). Azimuthal 
coverage for CLAS is limited by the magnet's six coils and 
is approximatively 90% at large polar angles and 50% at 
forward angles. 

The data were taken with an electron beam having an 
energy of 5.754 GeV impinging on an unpolarized 5-cm- 
long liquid- hydrogen target. The integrated luminosity of 
this data set was 28.5 fb~^. The data were taken from Oc- 
tober 2001 to January 2002. The kinematic domain of the 
selected sample corresponds approximately to Q^ from 1.5 
to 5.5 GeV^. We analyzed data with W, the 7* —p center- 
of-mass energy, greater than 1.8 GeV, which corresponds 
to a range of xb approximatively from 0.15 to 0.7. Here 

Xb is the standard Bjorken variable equal to w:i_^:i , qs 
with TUp the mass of the proton. 

The p° decays into two pions (tt+tt"), with a branching 
ratio of 100% [23]. To select the channel ep -^ e'pp^, we 
based our analysis on the identification of the scattered 
electron, the recoil proton and the positive decay pion 
(because of the polarity of the magnetic field, negative 
pions are bent toward the beam pipe and in general escape 
the acceptance of CLAS); we then used the missing mass 
ep —> e'pTr~^X for the identification of the ep —> e'pTr^Tr~ 
final state. 

Once this final state is identified and its yield normal- 
ized, the reduced 7*p -^ pp'^ cross section is extracted 
by fitting in a model-dependent way the (tt+tt^) invari- 
ant mass using a parametrized p'^ shape, which will be 
described later. The longitudinal and transverse cross sec- 
tions are then further extracted by analyzing the decay 
pion angular distribution in the p" center-of-mass frame. 
We detail all these steps in the following sections. 



2.1 Particle identification 

The electron is identified as a negative track, determined 
from the DC, having produced a signal in the CC and 
the EC. Pions, potentially misidentified as electrons, were 
rejected by cutting on the CC amplitude (> 2 photoelec- 
trons), imposing a minimum energy deposition in the EC 
(60 MeV) and correlating the measurements of the mo- 
mentum from the DC and of the energy from the EC. In 
order to minimize radiative corrections and residual pion 
contamination, a further cut E' > 0.8 GeV was also ap- 
plied, where E' is the scattered electron energy. Finally, 
vertex and geometric fiducial cuts, which select only re- 
gions of well understood acceptance, were included. 

The efficiencies of the CC and EC cuts, respectively 
Tjcc and r]Ec, were determined from data samples, select- 
ing unambiguous electrons with very tight CC or EC cuts. 
The CC-cut efficiencies range from 86 to 99% and the EC- 
cut efficiencies from 90 to 95%, depending on the electron 
kinematics. The efficiencies of the geometric fiducial cuts 



were derived from CLAS GEANT-based Monte-Carlo sim- 
ulations. 

Pions and protons arc identified by the correlation be- 
tween the momentum measured by the DC and the ve- 
locity measured by the SC. This identification procedure 
is unambiguous for particles with momenta up to 2 GeV. 
Particles with momenta higher than 2 GeV were there- 
fore discarded. The efhciencies of the cuts imposed for 
this momentum- velocity correlation and of the geometric 
fiducial cuts were determined from CLAS GEANT-based 
Monte-Carlo simulations. 

Once the electron, the proton and the positive pion 
are identified, the ep -^ e'pTT^'!r~ final state is identi- 
fied through the missing mass technique. Fig. [2] shows 
the square of the missing mass for the system e'p-K^ (i.e. 
M'^[e'pTT'^X]) as a function of the missing mass for the 
system ep (i.e. Mx[e'pX]). One distinguishes the p" and 
w loci quite clearly. A cut on the M'^[e'pTT^X] variable 
is required in order to separate the p'^ and the associated 
TT+TT" continuum from the u and the three-pion contin- 
uum background. The optimum value of this cut: 
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Fig. 3. Missing mass MJ.[e'p-n-+X] for W > 1.8 GeV and 
E' > 0.8 GeV. The red lines show the cut used (see eq. [ij to 
select the e' pn^ tt^ final state. 



mesons. These will be even more evident when we look at 
the differential spectra later on. 



0.05 < Mj.[e'pn+X] < 0.08 GeV^ 



(1) 



was determined from a study whereby we estimated the 
number of p° events, from fits to the M'^[e'pX] distribu- 
tion, as a function of the cut values. The cuts were chosen 
in the region where the number of p'^ events began to 
vary only very weakly with these cut values. The simula- 
tion used to calculate acceptances reproduces the features 
of fig. [3l The position of this cut relative to Al'^[e' pn'^ X] 
is shown in fig. [3] 





1.2 1.4 
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Fig. 4. Missing mass Mx[e'pX] for -0.05 < 
M%[e'p-K+] < 0.08 GeV^ W > 1.8 GeV and E' > 0.8 GeV. 
The p°(770), as weU as the /o(980) and /2(1270) resonances 
which can be distinguished, sits on top of a background of 
non-resonant two-pion continuum. 



Fig. 2. Squared missing mass Mx[e'piY'^X] vs Mx[e'pX] for 
W > 1.8 GeV and E' > 0.8 GeV. 



The missing mass distribution for the system ep, ob- 
tained after this cut, is shown in fig. [H The p" peak is 

-pO 



very broad : F^o « 150 MeV from ref. [23]. and sits on 
top of a background of a non-resonant two pion contin- 
uum, which originates from other processes leading to the 
e'pTr'^Tr~ final state, such as ep — > e'A^'^n'^ — > e'pT:'^n~ . 
In fig. m one can additionnally distinguish two bumps 
at masses around 950 MeV and 1250 MeV correspond- 
ing respectively to the scalar /o(980) and tensor /2(1270) 



2.2 Acceptance calculation 

The acceptance of the CLAS detector for the e'pn'^X 
process has been determined with the standard GEANT- 
based code developed for CLAS. Our event generator [53] 
contains the three main channels leading to the e'p'iT~^'K~ 



'A++ 



final state: ep -^ e'pp" ^^tt'tt ,ep— s-c'tt Zi' ' ^-s- pn ' 
and the non-resonant (phase space) ep — > e'pn'^Tr^ . This 
event generator is based on tables of total and differen- 
tial cross sections of double pion photoproduction data 



that have been extrapolated to electroproduction. This 
has been done by muhiplying these tables by a virtual 
photon flux factor and a dipole form factor in order to 
obtain a relatively realistic Q^ dependence of the cross 
section. We also have tuned the relative weight of all the 
aforementioned channels in order to reproduce the main 
kinematical distributions of our experimental data. 

Eight independent kinematical variables are necessary 
to describe a reaction with four particles in the final state. 
However, in unpolarized electroproduction, the cross sec- 
tion does not depend on the azimuthal angle of the scat- 
tered electron. The following seven variables are then cho- 
sen: Q^, ccb, i, M^+.;r-, ^1 cos{9hs) and (J)hs- Here Q^ 
and Xb are respectively the absolute value of the squared 
electron four-momentum transfer and the Bjorken vari- 
able, which describe the kinematics of the virtual pho- 
ton 7*. At some stages, W, the 7* — p center-of-mass en- 
ergy will equivalently be used. Then t is the squared four- 
momentum transferred to the p", <? is the azimuthal angle 
between the electron scattering plane and the hadronit 
production plane, and M^+^- is the invariant mass of the 
TT+TT^ system. Finally, cos{9hs) and (J)hs describe the de- 
cay of the p'^ into two pions and are respectively the pola: 
and azimuthal angles of the tt"*" in the so-called Helicit} 
Frame (HS) where the p° is at rest and the z-axis is giver 
by the p'^ direction in the 7* — p center-of-mass system 
All these variables are illustrated in fig. El 
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# of bins 


Width 


Q' 


GcV^ 


1.60 - 3.10 
3.10 - 5.60 
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Xb 
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0.16 -0.7 
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0.06 


-t 


GeV^ 


0.10 - 1.90 
1.90 - 4.30 
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0.30 
0.80 


<P 


deg. 


0.00 - 360.00 
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40.00 


cos(0^f) 
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-1.00 - 1.00 


8 


0.25 


J, -Hi' 


deg. 


0.00 - 360.00 


8 


45.00 


MMcp 


GeV 


0.22 - 1.87 


15 


0.22 



Table 1. Binning in the 7 independent variables for the ac- 
ceptance table. 




(p at rest) 



Fig. 5. Reference frames and relevant variables for the descrip- 
tion of the ep -^ e'pp° ^^ tt^tt" reaction. 



The procedure we have followed has been to calculate 
an acceptance for each of the T-dimensional bins. In the 
limit of small bin-size and unlimited statistics, this proce- 
dure is independent of the model used to generate events. 
Our event generator has nonetheless been tuned to repro- 
duce the experimental data. The binning in the 7 indepen- 
dent variables is defined in table [T] and its 2-dimensional 
(Q^, Xb) projection is shown in fig.[6l 

More than 200 million Monte-Carlo (MC) events were 
generated using CLAS GEANT to calculate the accep- 



Fig. 6. Binning in {Q ,xb) for our experimental data (with 
W > 1.8 GeV and E' > 0.8 GeV). 



tance in each of the individual 7-diincnsional bins. Each 
"real data" event was then weighted by the ratio of the 
number of MC generated to accepted events for each 7- 
dimensional bin. The acceptances are at most a few per- 
cent. Bins that have a very small acceptance (< 0.16%) 
have a very high weight, which produces an unphysically 
high and narrow peak in the weighted event distributions, 
and have to be cut away. The efficiency of this cut is evalu- 
ated by MC computation of the ratio of weighted accepted 
events to generated events mapped onto 1-dimensional dis- 
tributions. This correction factor ry^, is therefore model de- 
pendent since it is 1-dimensional and thus integrated over 
the remaining variables. It is on average 15%. 

Radiative corrections were part of our event generator 
and were calculated according to ref. [25] . The MC accep- 
tance calculation presented above therefore took into ac- 
count the effects of the radiation of hard photons and the 
corresponding losses due to the application of the cut of 
eq. [T] The contribution of soft photons and virtual correc- 
tions were determined by turning on and off the radiative 
effects in our event generator, defining an Frad factor for 
each (Q^, xb) bin for the integrated cross sections, or for 
each (Q^, xb, X) bin for the differential cross sections, X 
being one of t, (p, cos{9hs) or (f>HS- 



2.3 7*P -^ pTT+TT total cross section 

The total reduced cross section for the ep 
reaction can then be obtained from: 



to fully match each other. We will come back to this issue 
in section [275l when we are comparing the p° cross sections. 
e'pTT^TT^ The next step is to extract the j*p — > pp^ cross section 

from the j*p — > ptt+tt^ cross section, which requires a 
dedicated fitting procedure. 



^7*p--*p7r+7r IW :^Bi^) 



1 a (j e.p — >e'p7r+7r 



rv(Q ,XB,E) dQ dxB 2.4 Fitting procedure for the "f*p^pp^ cross section 



with: 



(T, 



ep — ye' pTT^TT 

dQ^dxB 



Lint AQ"^ AxB riccnECnw' 



(2) 



(3) 



where 



— nyj{Q^ ,xb) is the weighted number of ep —>■ e'pir^n^ 
events in a given bin {Q^, xb), 

— Lint is the effective integrated luminosity (that takes 
into account the correction for the data acquisition 
dead time), 

— AQ^ and Axb are the corresponding bin widths (see 
tabic [1]); for bins not completely filled, because of W 
or E' cuts on the electron for instance (see fig.O, the 
phase space AQ^Axb includes a surface correction and 
the Q^ and xb central values are modified accordingly. 

~ Frad is the correction factor due to the radiative effects 
(see section [2?2|) . 

— rjcc is the CC-cut efficiency (see section [^H]) , 

— rjEc is the EC-cut efficiency (see section [Ol) . 

— 77„ is the efficiency of the cut on the weight in the ac- 
ceptance calculation (see section [ 



We adopted the Hand convention [26] for the definition 
of the virtual photon flux Fy: 



rv{Q'',XB,E)^ 



1- Xb 1 



Stt to2£;2 



with 






(4) 
(5) 



and a « t4= the standard electromagnetic coupling con- 



stant. 



Fig.[7]shows the total reduced cross section 0-.^.^^^^+^- 
as a function of Q^ for constant W bins compared with 
the world's data [TSlfTOPTj- 

Relatively good agreement between the various exper- 
iments can be seen. It is important to realize that what 
is plotted is the unseparated cross section, i.e. a linear 
combination of the transverse (ctt) and longitudinal {(Tl) 
cross sections : a = ar + cctl. This means that, due to e 
(eq. [5]) , there is a dependence on the beam energy in this 
observable. Since the CORNELL data have been taken 
with an 11.5 GeV electron beam energy [TH], the DESY 
data with a 7.2 GeV electron beam energy [18] and the 
previous CLAS data with a 4.2 GeV beam |27], the data 
sets, although at approximatively equivalent Q^ and W 
values, are not directly comparable and are not expected 



Fig. [8] shows the acceptance- weighted Af^r+Tr- spectra for 
all our {Q^,xb) bins. The p'^ peak (along with the /o(980) 
and /2(1270) peaks clearly visible in some {Q^,xb) bins) 
sits on top of a tt+tt^ continuum background (see also 
fig. H] where all data have been integrated). 

This background can be decomposed, in a first approxi- 
mation, into the non-resonant ep — > e'pir^TT^ phase space 
and the exclusive electroproduction of a pion and a nu- 
cleon resonance, the latter decaying into a pion and a nu- 
cleon, such as ep — > e'n^A'^'^ ^-s- pn^ . Evidence for this 
can be seen in figs. [51 and [TUl which show for all our {Q^, 
Xb) bins the acceptance-corrected p7r~ andp7r+ invariant- 
mass spectra where structures are clearly seen. Most of 
these nucleon resonances (TV*) are rather well known, the 
most prominent being the Z\°'++(1232), the Di3(1520) 
and the Fi5(1680). However, their production amplitudes 
with an associated pion (i.e. ep —^ e'nN* '-^ p7r) are 
mostly unknown. 

At low energies {W <1.8 GeV) where very few N* 
can be produced, a phenomenological model has been de- 
veloped [^51l?7] based on an effective Lagrangian where a 
few N* 's are superposed along with the production of the 
p^. Such a model could be a strong constraint and guide 
to extract the p" cross section from all the other mecha- 
nisms. However, at the present higher energies, numerous 
new higher mass A^*'s appear as shown by the spectra of 
figs. [9l and [TOl For theoretical calculations, interference ef- 
fects between all these channels are virtually impossible to 
control and drastically complicate the analysis. Therefore 
this approach cannot be pursued in our case. 

At present, it is unrealistic to describe simultaneously 
the TT+TT", p'K~ and p7r+ invariant mass spectra over our 
entire phase space because there are too many structures 
varying independently with (Q^, xb) in each invariant 
mass distribution. 

Therefore, since a complete description of the di-pion 
mass spectra is not available, we have adopted an em- 
pirical description of the data using non-interfering con- 
tributions that, together with a model for the p° shape, 
reproduce the tt+tt^ invariant mass spectrum. 

The p° peak is broad and the strength of the non- 
resonant TT+TT^ background under it is quite significant, 
and, even more importantly, its nature is unknown. There- 
fore we must carry out a non-trivial and model-dependent 
fitting procedure in order to extract the 7*p — > pp^ cross 
section. 

The procedure we have followed consisted of fitting 
only the M^-^^- distributions for each {Q^, xb) bin in 
the case of the integrated cross sections and for each (Q^ ^ 
Xb, X) bin, where X can be t, <?, cos.{6hs) or (/)ifs, in 



X10 

k 

+ 

Q. 

t 
Q. 

'^ 1 



+ 

Q. 

T 

Q. 

'^ 1 

D 



1.80 <W (GeV) < 2.00 


2.00 <W (GeV) < 2.20 


2.20 <W (GeV) < 2.40 


* 


i-n-H f 

* 


T * 




2.40 < W (GeV) < 2.60 


2.60 <W (GeV) < 2.80 




: o 

. D^ 

, . , 1 . , , . 1 . , . , 1 . , , , 1 . , . , 1 . , . , 1 . , . 


■ 1 ■ 1 1 ■ 1 1 1 1 1 ■ 1 1 1 ■ 1 ■ 1 1 1 1 ■ 1 ■ 1 ■ 1 1 1 1 1 










n CORNELL 

o DESY 

A CLAS (4.2 GeV) 

• CLAS (5.754 GeV) 



2 3 4 5 6 
Q^ (GeV^) 



2 3 4 5 6 
Q^ (GeV^) 
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the case of the differential cross sections. It was possi- 
ble to fit the Mjr+T- spectra with five contributions: three 
Breit-Wigner shapes to describe the three evident mesonic 
TT+TT" resonant structures of the p°{770), /o(980) and /2(1270), 
where the masses in MeV indicated in parentheses are 
the values given by the Particle Data Group (PDG) p3] . 
and two smoothed histograms that are the A/^r+Tr- pro- 
jections of the reactions ep — *■ e'7r^Z\++ ^-> pn^ and of 
the non-resonant continuum ep —^ e'pn^n^. These two 
latter spectra are calculated by our aforementioned event 
generator [53] . We now detail these five contributions and 
explain why they are necessary (and sufficient). 

We first discuss the contribution of the p°(770) and 
the way to model it. It is well known that simple sym- 
metric Breit-Wigner line shapes which are, to first order, 
used to describe resonances, are too naive to reproduce 
the p shape, because of, among other aspects, interference 
effects with the non-resonant tt+tt" continuum. Several 
methods can be found in the literature for treating the p*^ 
shape (see for instance ref. [5S] for such a discussion) . The 
procedure we adopted was the following: 



— Introduction of an energy-dependent width in order 
to take into account that the p^ is an unstable spin- 
1 particle that decays into two spin-0 particles; it is 
also called a p-wave Breit-Wigner 0^. This modified 



Breit-Wigner reads: 

BWp{M,+^- 



Mpr{M^+^-) 



(A/2-M2+^_)2-fAf2r2(M.+.-) 



with the energy-dependent width: 

21 + 1 






M„ 



(6) 



(7) 



where I = 1 for a p-wave Breit-Wigner, q is the momen- 
tum of the decay pion in the p'^ center-of-mass frame 
and qp is equal to q for A/^+j^- = A/^: 



9 = 



M' 



4Af2 



9p = 



Af 2 - 4Af 2 



(8) 



Ross and Stodolsky [3T| and Soding 32J have shown 
that the interferences between the broad p^ peak and 
the important non-resonant tt+tt" contribution under- 
neath leads to a skewing of the Breit-Wigner. Accord- 
ing to Ross-Stodolsky, one way to take account of this 
effect is to introduce a correction term that consists of 
multiplying the Breit-Wigner formula by an empirical 
factor that shifts the centroid of the A/^+^- distribu- 
tion: 



BW;''-{M,+^-)^BWp{M,+^- 



M„ 



M^+^- 



(9) 
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at Mx ~ 0.770, 0.980, 1.275 GeV corresponding to the three well-known resonant states in the (tt^tt^) system. 



where Uskew is the "skewing" parameter. Although 
Ross-Stodolsky have predicted the value of Uskew to 
be 4, it is often a parameter that is fitted to the data 
since so little is known concerning the interference be- 
tween the p" signal and the tt+tt" continuum. 

As is evident in fig. [3 in addition to the /9"(770), there 
are two well-known resonant structures in the tt+tt^ sys- 
tem: the /o(980) and /2(1270). Due to the large widths 
of these mesonic resonances (40 to 100 MeV [23] for the 
/o(980) and ^i 180 MeV l23j for the /2(1270)), it is clearly 
necessary to include them in our M^+^- fit because their 
contribution can extend into the p^ region, which is itself 
also broad. We also used the formulas of eqs.[6][9]for these 
two other mesonic resonant states /o(980) and /2(1270) 
with appropriate parameters Mf^ , Ff^ , Mj^ , Ff^ and took 
into account their I — and I — 2 nature. 

In principle, the only free parameter to vary in eq. [9] 
should be ngkew However, we have also allowed the cen- 
tral masses and widths of the three mesons to vary in a 
very limited range of at most 20 MeV from their nominal 
values (see table [2]). The motivation for this is that, be- 
sides the largely unknown interference effects between the 
meson and the tt+tt" continuum, several other effects can 



shift or distort the meson shapes: radiative corrections, 
binning, acceptance corrections, imprecise values for the 
central masses and widths of some of these mesons, etc. 

Finally, besides the p°, /o(980) and /2(1270) mesons, 
the two other contributions entering our fit are: the M^+^- 
projections of the reactions ep —> e''K~ A^~^ ^-> pTr"*" and 
the non-resonant continuum (phase space) ep — )■ e'pvr+Tr^. 
The shapes of these distributions are given by our event 
generator. In particular, the Z\++ has the shape of a stan- 
dard Breit-Wigner in the pTr"*" distribution with a centroid 
at 1.232 GeV and a width of 111 MeV [23]. As is obvious 
from fig. [9l the A^ contribution can be neglected. 

In principle, of course, other processes contribute to 
the M^+T^- continuum, for instance all of the ep —^ eTr~ N* ^ 
p7r+ reactions, as already mentioned. As a test, we mod- 
eled the p'K~ and pTr"*" invariant mass distributions of figs. [9] 
and [TUl bv adding, at the cross section level, several Breit- 
Wigners matching the structures seen in these figures and 
identifying them with known N* masses (and widths) that 
can be found in the PDG. Like for the Z\"'""'", we introduced 
their contribution into our fit of the M^+^r- distribution 
using our event generator. The conclusion we reached was 
twofold. Firstly, this procedure introduced a large number 
of additional free parameters: for each of the extra iV*'s, 
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parameter 


PDG value 


min. value 


max. value 


p" mass Mp (MeV) 


^ 770 


750 


790 


p" width Fp (MeV) 


Ri 150 


140 


170 


/o mass Mfg (MeV) 


^980 


970 


990 


/o width T/o (MeV) 


40-100 


40 


120 


/a mass Mf^ (MeV) 


^ 1275 


1260 


1280 


/a width Ff^ (MeV) 


^ 185 


170 


200 



Table 2. Range of variations permitted for the parameters to be fitted in formula [9] 



two parameters for the central mass and width to vary in 
the approximate ranges given by the PDG and one more 
for the weight/normaUzation. Secondly, we found that the 
-^TT+TT- projected shape of these high mass iV*'s was very 
similar to the phase-space Mj^+^^ distribution. In other 
words, in a first approximation, the phase-space contri- 
bution can reflect and absorb the high mass A^*'s. How- 
ever, the Zi++ contribution to the Mj^+^- distribution was 
found to be sufficiently different from the phase space dis- 
tribution to be kept as an individual contribution. 

To summarize, each {Q^, xb) bin of fig. [H] was fit with 
the following formula: 



dN 



BWp{M^+^-) 

+BWf„ (Af,+,- ) + BWf, (M,+,- ) 
+Af^++,- (M,+,- ) + Mp,+,- (A'4+,- ) 

(10) 



It involves 14 parameters that are: 

— 1) weight (normalization), 2) central mass, 3) width 
and 4) n^kew oi p°; 

— 5) weight (normalization), 6) central mass, 7) width 
and 8) Uskew of /o; 
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9) weight (normalization), 10) central mass, 11) width 
and 12) Uskew of /2; 

13) weight (normalization) of Mtt+tt- projection of the 
ep —^ e'7r~Z\ "'""'" ^-> pvr^ process; and 

14) weight (normalization) of ep — > e'pTT^n^ phase 
space. 



Fourteen parameters might appear a lot to fit only a 
1-dimcnsional distribution. However, on the one hand, six 
of these (the central mass and width of the /9°, /o(980) and 
/2(1270) mesons) are quite constrained and are allowed to 
vary in a very limited range. On the other hand this simply 
reflects the complexity and our lack of knowledge of the 
ep —^ e'pT:'^TT~ reaction, to which many unknown, inde- 
pendent though interfering processes contribute; namely: 
meson production ep -^ e'pM^ ^^ tt+tt", TV* production 
N*++ ^ nTT+, en -^ e'TT+N*° 



ep 



e TT 



pn 



resonant ep —^ e ptt+tt , etc. 

Fig. fTTJ shows the result of our fits to the M^+^- distri- 
butions, normalized in terms of the reduced cross sections 
of eq. [3] for all of our (Q^, xb) bins. In a few cases, the 
fits do not fully describe the data. For instance, for 0.46 
< Xb < 0.52, 3.10 < Xb < 3.60, the data tend to show 
a "structure" around M^+jr-=0.9 GeV, i.e. between the 



known p*^ and /o resonances, which cannot be reproduced 
by our fit formula of eq. [TOl We attribute this discrepancy 
to interference effects not taken into account by our sim- 
ple fit procedure. As discussed in more detail in the next 
subsections, a systematic uncertainty of 25% is assigned 
to this whole fit procedure which is meant to account, 
among other aspects, for the inadequacies in the model. 
On all the figures that are going to be presented from now 
on, unless explicitely stated otherwise, all the error bars 
associated to our data points will represent the quadratic 
sum of the statistical and systematic errors. 



2.5 Integrated p° cross section 

We use the p'^ strength (green line) extracted from the 
distributions shown in fig.[TT]to calculate the cross section. 
Fig. [T2l shows the resulting reduced cross section cr^'p^ppo 
compared with the world's data presented as a function 
of W for constant Q^ bins. Fig. [13] shows the reduced 
cross section a^*p^ppO compared with the world's data 
presented as a function of Q^ for constant W bins. 

With respect to the 7*p -^ ptt^tt^ cross section that 
we extracted in the previous section, there is an additional 
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7 P->PP 



,0 cross 



source of systematic uncertainty for the a. 
section that arises from the subtraction procedure de- 
scribed in the previous section. This contribution is quite 
difficult to evaluate. It is not so much the quality of the 
fit in fig. 111! that matters; we have varied the minimum 
and maximum limits imposed on the parameters in ta- 
bleland found that the results of the fits are very stable. 
The uncertainty arises more from the reliability and confi- 
dence we can assign to the modeling that we have adopted 
for the p'^, /o(980) and /2(1270) mesons with the skewed 
Breit-Wigners and for the non-resonant continuum tt+tt" 
distribution. We have tried several shapes for this latter 
continuum. As mentioned in the previous subsection, we 
introduced N* states other than the Z\++. Ultimately, we 
ended up finding the fits to be stable at the « 20 % level 
on average. Overall, we cannot take account of any inter- 
ference effects between the p'^ peak and the non-resonant 
7r+7r~ continuum. This uncertainty is of a theoretical na- 
ture, and in the absence of sufficient guidance at present, 
we have decided to assign a relatively conservative 25% 
systematic uncertainty to our extracted p" yields. We will 
find some relative justification for this estimation in the 
next section when we study the differential distributions, 
in particular those of t and cos 9 hs- 

Recently, a partial wave analysis of data of exclusive 
TT'^TT" photoproduction on the proton from CLAS, has 
been carried out [34l . This study showed that the p° cross 



sections resulting from this sophisticated method were 
consistent with those resulting from simple fits of the two- 
pion invariant mass as we have just described, to a level 
much lower than 25%. Although no such partial wave anal- 
ysis has been done to the present electroproduction data, 
this photoproduction comparison gives relative confidence 
that the 25% systematic uncertainty that we presently as- 
sign, is rather conservative. 

Coming back to fig. [121 we find that our data are in 
general agreement with the other world's data in regions 
of overlap. In the upper left plot of fig.[l2](1.60< Q^ <i.90 
GeV^), our CLAS (5.754 GeV) data seem to overestimate 
the CLAS (4.2 GeV) results, but this can certainly be 
attributed to a kinematic effect due to the different beam 
energies of the two data sets. Indeed, we are comparing 
the total reduced cross sections : a — ax + ^ctl- However, 
at M^=2.1 GeV and Q^=l.7 GeV^, e=0.53 for a 4.2 GeV 
beam energy but e=0.77 for a 5.754 GeV beam energy. 
This can readily explain the lower CLAS (4.2 GeV) data 
with respect to the CLAS (5.754 GeV) data. 

On this general account, we could have expected that 
the Cornell data stand to some extent above the CLAS 
(5.754 GeV) data since they have been obtained with an 
11.5 GeV beam energy. This is not the case which might 
indicate a slight incompatibility between the Cornell and 
CLAS data. This point, as well as the compatibility of 
the CLAS (4.2 GeV) and CLAS (5.754 GeV) data, wiU be 
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confirmed in section [3?al wliere we compare tlie separated 
longitudinal and transverse cross sections for which this 
beam energy kinematical effect is removed. 



2.6 Differential p" cross sections 

After having obtained the total p° cross section, we now 
extract the differential cross sections in t, <?, cos{9hs) and 

Since the data are now binned in an additional vari- 
able, each bin has fewer statistics, not only for the real 
data but also for the MC data that are necessary to cal- 
culate the acceptance correction. Bins for which r]yj is less 
than 0.6 were rejected, where, we recall, rj^ is the correc- 
tion factor in the acceptance calculation that was intro- 
duced in section [2T2l This explains why some holes occur 
at several instances, in particular in the <P and cos{9hs) 
distributions. 

We start by extracting the da/dt cross section. Defin- 
ing i' as t — to, where to is the maximum t value kine- 
matically allowed for a given (Q^, xb) bin, we divided the 
data into 6 bins for < — i' < 1.5 GeV^ and 3 bins for 
1.5 < -t' < 3.9 GeV^ For each of the (t, Q^, xb) bins, we 
extracted the p^ signal from the (7r+,7r^) invariant mass 
spectra using the fitting procedure previously described. 

Fig. [Ml shows da/dt for all our {Q'^,xb) bins as a func- 
tion of t. The general feature of these distributions is that 
they are of a diffractive type, i.e. proportional to e''*. The 



values of the slope b are between and 3 GeV~^. They 
are plotted as a function of W in fig. [15] along with the 
world's data. For the sake of clarity, only the world's data 
for Q2 > 15 GeV^ are displayed. For Q^ < 15 QeV^, 
the data show the same trend but with more dispersion. 
The data exhibit a rise with W until they reach a plateau 
around W= 6 GeV at a 6 value of « 7 GeV'^. The high- 
energy experiments (HI and ZEUS) have shown that this 
saturating value tends to decrease with Q^, which is il- 
lustrated by the HI points in fig. [TSj that correspond to 
different Q^ values. 

By integrating the da/dt cross section, we are able to 
recover at the « 20% level the integrated cross sections 
that were presented in section 12.51 The agreement is not 
perfect since for the integrated cross section one fits a sin- 
gle full statistics M^+^- spectrum, whereas for the differ- 
ential cross section, one fits several lower statistics M^+^- 
spectra, that are then summed. This relatively good agree- 
ment serves, among other arguments, to justify the 25% 
systematic uncertainties that we have applied in the non- 
resonant 7r+7r~ background subtraction procedure. 

We note that the integrated cross sections that we have 
presented so far (and which will be presented in the re- 
mainder of this article) have been summed over only the 
domain where we had data and acceptance. We have not 
extrapolated our cross sections beyond the t domain ac- 
cessed in this experiment, which we deem unsafe and very 
model-dependent. Fig. [T3| indicates that this might under- 
estimate some integrated cross sections for a (very limited) 
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number of {Q'^,xb) bins at large xb, where the t depen- 
dence appears rather flat. 

We proceeded in the same way to extract da/d'P. All of 
our {Q^,xb) bins are shown in fig. [El Several of the bins 
near ip—lSO" are empty or have large error bars because 
of very low acceptance in CLAS in this region. 

These distributions were fitted with the expected <? 
dependence for single meson electroproduction: 



da 
d^ 



2tx 



( ot + foi^ 



+ecos2<? ott + A/2e(l + e)cos^ otl) 



(11) 

from which we could extract the interference terms gtt 
and GxL ■ The curves in fig. [12] show the corresponding 
fits, and ot + eci,, utt and otl are displayed in fig. [T71 
If helicity is conserved in the s channel (SCHC), the in- 
terference terms ott and aj^j^ would vanish. Most of our 
extracted values are consistent with within (large) er- 
ror bars, although one clearly cannot make strong claims 
about SCHC at this point. 

Turning to the pion decay angles of the p°, dus and 
4>HS, they are expected to follow the general and model 
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independent distribution |40| 
VF(^, cos 6*^5, </?ffs) = 

- V2Rer°o sin 2Bbs cos ipus - r^x-x sin^ Bus cos lipus 

-e cos 2<?(rJi sin^ dus + ''oo cos^ Ors 

-\/2Rerio sin 2dns cos kpns ~ A-\ sin^ Qhs cos Ikpus) 

— e sin 2<?( v2ImrjQ sin IQhs sin ipHs 

^-\mr\_^ sin^ Qus sin 2(pHs) 



V2e(l + e)cos<?(r^i sin^ Bus + r^o cos^ 61^5 
- V2Rerio sin 26*^5 cos if as - t\_\ sin^ Ors cos 2kpHs) 
-^26(1 + e)sin ^(V2ImriQ sin 2Qb_s sin if us 
-Imr^_i sin^ Brs sin "l^Hs) \ , (12) 



where: 



„04 



P: 



eRpjj 



1^ 

P% 



eR 



1 + eR 



1,2 



15 



> 

0) IO- 
CS '"i^ 



6- 



• CLAS (5.754 GeV) 

a CORNELL (1979,1981) 

HERMES 

* NMC 



A FermiLab(1977) 
o E665 
*c H1 
• ZEUS 




J I I I L_ 



J I L_ 



± 



10 



W(Ge 



i% 



Fig. 15. Slope b of da/dt as a function of W. Data of Cornell [T9]|35] , HERMES [20], NMC [36], Fermilab (1979) [37], E665 [33], 
HI [38] and ZEUS [39] are shown for comparison. 



= Vi? 



1 + eR 



a = 5,6 



(13) 



with Rp equal to the ratio aj^jaT- 

The parameters pf- are bilinear combinations of the he- 
licity amplitudes that describe the 7*p -^ p'^p transition. 
They come from a decomposition of the 3x3 spin den- 
sity matrix of the p^ in a basis of 9 hermitian matrices. 
The superscript a refers to the virtual photon polariza- 
tion: a ~ 0-2 for transverse photons, a = 4 for longitu- 
dinal photons, and a = 5-6 for the interference between 
L and T terms. The subscript refers to the vector meson 
helicity: i,j — refers to a longitudinal polarization state 
and I, j = —1, 1 to a transverse polarization state. For ex- 
ample, Pqq is related to the probability of the transition 
between a transverse photon (a = 0) and a longitudinal 
meson (i,j — 0) and p^i is an interference term between 
meson helicities and 1 (i — 0,j ~ 1) produced by a 
transverse photon (a — 0). 

If SCHC applies, then by definition, pgp = and Pqq = 
1. Then eg. [T3l leads to a direct relation between the mea- 
sured Tqq and the ratio Rp — ^^. In that case, the longi- 
tudinal and transverse cross sections, cr^ and ax, may be 
extracted from the cos 0hs distribution, without relying 
on a delicate Rosenbluth separation. 

SCHC can be tested by studying the integrated distri- 
butions W{(j>Hs) and W{'P) (where 'P — cfiHS ~ ^) over 
the other decay angles. Integrating W{(P, cos9hs, 4>hs) of 



eq. [12] over cos Ohs and <? yields 
1 



W{(j) 



HS 



2tt 



[l-2r°iicos20ffs], 



(14) 



which isolates rjlj^, a density matrix element violating 
SCHC. Integr a.tmg W{<1^, cos 9 HS,4'Hs) over cos 0hs yields: 



Wi'P) = -!- [l + 2erl_^ cos2tZ>'] 



(15) 



Another consequence of SCHC is that the W{(I)hs) 
distribution should be constant and the W{'1') distribution 
should vary as cos 2'!' if rl_i is not zero. 

We extracted the da/d(j)HS and da/dP cross sections 
in the same way as previously mentioned, i.e. by fitting 
the Mx [e'pX] spectra and extracting the p^ yield for each 
{Q^,xb,(I)hs) and {Q'^,xb,^) bin, respectively. Fig.[T5lshows 
the extracted values of r^lj^ and rl_i, obtained by fitting 
da/d(i)HS and da / dP with the functions of eqs.fMlandlTSl 
respectively. Basically all the SCHC violating matrix el- 
ements r^'^i are compatible with (though within large 
uncertainties) , which gives some relative confidence in the 
validity of SCHC. In addition, r\_^ is also found to be 
compatible with for all kinematics, although this is not 
a necessary requirement for SCHC. This indicates that our 
P distributions are basically flat. 

We finally extract the Tqq matrix element from the 
cosOhs distributions, which result from the integration of 
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0.70 



W{(1', cos 9hs, 4>hs) of eq. [T^lover (f>Hs and <P: 



Wicosens) = ^ [(1 - ^o"o) + (3r[l^ - 1) cos^ 



^HS 



(16) 



As an example, fig. [19] shows a cos 9hs distribution, 
before and after the non-resonant 7r"'"7r~ background sub- 
traction, for one of our typical {Q'^,xb) bins. We note that 
the unsubtracted distribution is highly asymmetrical in 
cos Ohs ■ This is mainly due to the presence of events from 



the ep 



'/i++7r reaction, whose phase space is maxi- 



mum around cos 6hs = 1 • 

Fig. HOI shows the da / d cos{9 h s) cross sections for all 
bins in {Q^,xb)- Even after the non-resonant 7r"^7r~ back- 
ground subtraction procedure, some of the aforementioned 
asymmetry in the cos{9hs) distribution remains at the 
« 25% level. We attribute this to interference effects be- 
tween the p'^ channel and its background (mostly, ep —^ 
e'n^ A*^^ ^^ p7r+ and non-resonant ep — > e'pTT^n^ as 
already discussed), which obviously cannot be taken into 
account when subtracting the different channels at the 
cross-section level as we do. This ~ 25% cos 6hs asymme- 
try is a further confirmation of the systematic uncertainty 
associated to the extraction of the p^ signal. Fig. [T51 shows 
then the resulting rgp values. 
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Fig. 19. An example of a (acceptance corrected) cosOhs 
distribution before (black points) and after (red) the non- 
resonant -k'^-k~ background subtraction (0.58 < xb < 0.64 
and 4.10 < Q^ < 4.60 GeV^ bin). In this example the error 
bars are purely statistical and no systematic uncertainty has 
been added. The asymmetry in the red data points between 
cos&_r/s=-l and cos6'hs=1 is attributed to some remaining 
non-resonant tt^tt" background which could not be subtracted 
by our fitting procedure, estimated to lead to a 25% systematic 
uncertainty (see section [274)) . 
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2.7 Longitudinal/transverse cross section separation 



Following the relative verification of the presence of SCHC 
in the previous discussion, the ratio Rp can be determined 
from: 



The separated longitudinal and transverse cross sec- 
tions, (Tl and ctt, are then calculated as: 



(Tt = 



(1 + eRp 



(^L = 



{l + ^RpY 



(19) 



Ro 



(7L 

(Jt 



1 rl 



,04 



£ ^ '00 



(17) 



Although we cannot claim that our data give strong 
evidence for SCHC (nor for its violation), it should be 
noted that ref. ^5T] mentions that eq. [T7] is relatively ro- 
bust to violations of SCHC. Fig. [^T] shows Rp for all our 
{Q^,xb) bins, assuming SCHC, which is an assumption 
that we will keep for the remainder of this analysis. 

We fit our 27 points to a linear function: 



Rp = a + bQ^, 



(18) 



which yields: a = 0.281 ±0.549 and b = 0.439 ±0.203. The 
uncertainties on a and b are relatively large but highly 
correlated. The normalized correlation coefficent is 0.966. 
Fig. [221 shows that the band corresponding to our fit is in 
good agreement with the world's data. 



We will display these resulting cross sections in sec- 
tion [3l where they will be compared to the theoretical 
models. 



2.8 Systematic uncertainties 

There are several sources of systematic uncertainties. The 
main one stems undeniably from the fit to extract the p^ 
cross section and subtract the non- resonant TT'^n~ contin- 
uum, as discussed in section [2^ and estimated to be 25%. 
We recall that this estimation arises from several analyses: 

— Fitting the M^^+j^- distributions of fig. [TT] and chang- 
ing the shapes of the various inputs and contributions. 
Other systematic studies on these fits were carried 
out: removing a few data points on the edges of the 
M^+^- spectra to study edge effects, smoothing the 
histograms to take into account potential statistical 
fluctuations, varying the ranges of the parameters to 
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be fitted (see Table [2]), etc. All in all, we found a sta- 
bility and robustness of our fits at the w 20% level. 
Integrating da/dt over t (section I2.6P and comparing 
it to the total cross sections f section [2.5|) . resulting in 
a « 20% level agreement. 

Observing an asymmetry at the « 25% level between 
the forward and backward angles in the cos Ohs distri- 
butions (section [ 



A second source of systematic uncertainty stems from 
the acceptance calculation which is largely model indepen- 
dent. We have carried out several tests to determine the 
associated uncertainty on our procedure. For instance, we 
have varied the binning of the 7-dimensional table (see ta- 
ble[T|). As an other test, we have also varied the input event 
generators: taking for instance only the ep -^ e'pT:^TT~ 
phase space channel or only the ep ~~* e'pp^ channel. Ul- 
timately, we estimated the stability of our acceptance to 
be at the 15% level. 

We correct the data for radiative effects. These were 
generated according to ref. |25j . The approximations used 
in this calculation may lead to systematic uncertainties 
that we estimate to be of the order of 4%. 

The determination of the CC efficiencies relies on the 
assumption that the distribution of photoelectrons for de- 
tected electrons is a generalized Poisson function and that 
the shape of this distribution above 4 photoelectrons is 
sufficient to determine the whole distribution and to ex- 
trapolate to 0. The maximum error we could make on the 
integral of the distribution between and 4 photoelectrons 
is about 25%. Since the CC inefhciencies are at most of 
the order of 6%, the corresponding systematic uncertain- 
ties on the cross sections is 1.5%. 

The EC efficiency determination relies on the assump- 
tion that a particle with a sufficiently high number of pho- 
toelectrons in the CC was unambiguously an electron. We 
tried applying several values for this cut to estimate the 
stability of the results and the maximum differences were 
found to be on the order of 2%. 

The target length is known to about ± 1 mm. The 
hydrogen density was kept fairly constant through tem- 
perature and pressure stabilization. The determination of 
the beam integrated charge also has a small systematic 
uncertainty. All this is summarized in table [H and leads 
to a normalization error applicable to the whole data set. 



Source of error 


Estimated uncertainty 


Fitting procedure 


25% 


CLAS acceptance 


15% 


Radiative corrections 


4% 


CC efficiency 


1.5% 


EC efficiency 


2% 


Target thickness 


2% 


Target density 


1% 


Beam integrated charge 


2% 



Table 3. Systematic uncertainties affecting the overall nor- 
malization. The quadratic sum of all theses errors results in a 
~ 30% systematic error bar. 



3 Theoretical interpretation 

3.1 The Regge "hadronic" approach 

The Regge approach consists of understanding exclusive 
p" electroproduction, above the resonance region and at 
forward angles where the cross section is the largest, in 
terms of exchanges of meson "trajectories" in the i-channel. 
Regge theory generalizes the notion of a i-channel single 
particle exchange to the notion of a family (i.e. trajectory) 
of particle exchanges. Indeed, mesons, and more broadly 
hadrons, appear in general in sequences made of rotational 
excitations. Mesons that have the same quantum numbers, 
except for spin, seem to align along linear "trajectories" 
a{t) = a{0) + a't that relate their squared mass —t to 
their spin a. This leads in the high energy limit to ampli- 
tudes proportional to s"^*\ where s — W^, and therefore 
total cross sections are proportional to s"(°)^^. 

In the following, we will use "JML" to refer to the lat- 
est version [42^ of the model developed by J.-M. Laget and 
collaborators [TTlfT2[|43] . The dominant amplitudes corre- 
spond to the i-channel exchange diagrams of fig.[23l Since 
vector mesons have the same quantum numbers as the 
photon, systems with quantum numbers of the vacuum 
can be exchanged. The corresponding trajectory is called 
the pomeron. Although the pomeron contribution fully ex- 
plains meson photo- ^44] and electroproduction [45], it 
represents only about one third of the cross section in the 
1*P ~^ PP^ channel for the energy range covered by our 
data. Here the bulk of the cross section comes from the 
exchange of the /2(1270) and a mesons. The exchange of 
the TT meson, which dominates lj production, contributes 
very little to the p° production channel. 

Amplitudes for the pomeron, tt and /2 meson exchange 
diagrams can be found in ref. [11] and for the a me- 
son exchange diagram in ref. [43j . In photoproduction, 
the only parameters of the model are the coupling con- 
stants at the vertices of the diagrams. They are taken 
from a comprehensive study of other independent pro- 
cesses. For instance, the quark-pomeron coupling constant 
is fixed by the analysis of pp scattering at high energy 
[W ^ 100 GeV). In electroproduction, a monopole form 
factor is introduced at the ^np vertex and a dipole from 
factor is used at the "fap vertex [TTj. In this latter refer- 
ence, a dependence on t is given to the cut-off mass that 
accounts for an increasing point-like behavior of the cou- 
pling of the photon with the meson when —t (and conse- 
quently the impact parameter) increases. By construction, 
both the Q^ and the t dependency of the 7pp (where p 
stands for the pomeron) and the 7/2/0 vertices are intrin- 
sically part of the corresponding amplitudes, and no other 
parameters are included. 

With this limited number of parameters, the JML model 
is able to successfully reproduce the main trends of the W 
and t dependences of the total and differential cross sec- 
tions for the reactions 7*p — > pip, 7*p -^ puj and j*p — > 
pp^ over the whole W range, i.e. from threshold up to 
HERA energies. In order to save computation, in our case, 
the pomeron exchange version of the model has been used 
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Fig. 23. Dominant f-channel exciiange diagrams for tlie reac- 
tions 7p -^ pp^ . 



Fig. 24. Tiie liandbag diagrams for (longitudinal) vector me- 
son production. Quark GPDs are accessed on the left and gluon 
GPDs are accessed on the right. 



instead of the two gluon exchange version. In the momen- 
tum transfer range of this work, i.e. —t<2 GeV^, the 
two models lead to almost identical results. Overall, as 
was mentioned before, the pomeron/two-gluon exchange 
contribution does not dominate the cross section in the 
energy range that is accessed in this study {W up to 2.5 
GeV). 



3.2 The GPD "partonic" approach 

The JML model was originally built for photoproduction, 
i.e. Q^ = 0, and was extended to electroproduction by 
introducing form factors to take into account the shorter 
distances probed by the virtual photon, inversely related 
to Q^. We now consider another approach, based on the 
formalism of Generalized Parton Distributions (GPDs), 
which is valid in the so-called Bjorken (or "Deep Inelas- 
tic") regime, i.e. Q^ji' —» oo with xb — 2mI^ finite. An 
important question is how low in Q^ this asymptotic for- 
malism can still be applied or extrapolated? 

Collins et al. [7] have shown that the dominant pro- 
cesses for exclusive meson electroproduction, in the Bjorken 
regime, are given by the so-called handbag diagrams rep- 
resented in fig. [M] 

The handbag diagrams are based on the notion of fac- 
torization in leading-order pQCD between a hard scatter- 
ing process, exactly calculable in pQCD and a nonpertur- 
bative nucleon structure part that is parametrized by the 
Generalized Parton Distributions introduced by Miiller et 
al. [1], Ji [5] and Radyushkin [6]. For the quark handbag 
diagram and for vector mesons, only the two unpolarized 
GPDs contribute. They are called, using Ji's notation, H 
and i?, and they depend upon three variables: x, ^ and t. 
We refer to the rich literature on GPDs (see refs. piQlfTO] 
for recent reviews) for the full definition of the formalism 
and of the variables. For the gluon handbag diagram, the 



corresponding H^ and E^ gluonic GPDs are usually ap- 
proximated by and reduced to the forward gluon density 
G{x). 

We recall that for mesons, factorization, which is an 
essential component of the handbag mechanism, is only 
valid for the longitudinal part of the cross section, as the 
L subscripts on the photon 7^ and on the meson p° in- 
dicate in fig. [211 This is one of the main motivations for 
separating the longitudinal and transverse parts of the 
cross sections in our data analysis. 

We stress that theoretical calculations of exclusive me- 
son production cross sections in the QCD factorization 
and GPD approaches are extremely challenging because 
one needs to address several issues at once: how to model 
the GPDs; how to treat the hard scattering process (the 
choice of an effective scale in as, the role of QCD cor- 
rections, etc.); how to consistently combine contributions 
from meson production in small-size and large-size con- 
figurations, etc. While in theory these are distinct issues 
that can be discussed separately, in practice they are very 
much related. Therefore, the choices and approximations 
one makes in the treatment of one will generally infiuence 
the conclusions one draws about the others. 

In the following, we will discuss the two particular 
GK [13] and VGG [M l fTS l fTeT?^ GPD-bascd calculations 
that provide quantitative results for the longitudinal ex- 
clusive p° cross section. Both groups have adopted the 
same approach. They parametrize the {x,S,) dependence 
of the H and E GPDs based on double distributions as 
proposed in ref. [46] , (the treatment of the t dependence 
being different). They correct the leading order amplitude 
with an intrinsic transverse momentum dependence, the 
so-called k± corrections or, more generally, the modified 
perturbative approach [47j . On this latter point, it is in- 
deed well-known that at high W the leading-twist calcu- 
lations overestimate the data and that the associated pre- 
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diction of ^^t evolving as 755- , at fixed xb , is not observed 
in the data (see for instance ref. [48j). 

The main difference between the two calculations lies 
in the fact that the GK group has treated the sum of the 
two handbags at the amplitude level, while the VGG group 
has treated it at the cross section level, and has therefore 
neglected the interference between the two handbag dia- 
grams. We will see this effect in the next section where we 
compare our data to the two particular GPD models we 
have just introduced. 



3.3 Comparison to data 

Fig. [25] shows our results for the total longitudinal cross 
section for exclusive p'^ electroproduction on the proton 
as a function of W, for different Q^ bins, along with the 
relevant world's data. 

The cross sections clearly exhibit two different behav- 
iors as a function of W. At low W the cross sections de- 
crease with W and then begin to rise slowly at M^ « 10 
GeV. 

In fig. [25] the results of the calculations of the JML, 
VGG and GK models are also shown. The JML model 
(dash-dotted line) reproduces fairly well the two general 
behaviors just mentioned. The drop of the cross section at 
low W is due to the i-channel a and /2 meson exchange 
diagrams (see fig. [23]). The intercept q;(0) of the /2 trajec- 
tory is w 0.5 and therefore the cross sections decrease with 
energy as -^ . The flattening of the cross section near W ~ 
10 GeV comes from the combined effect of the decreasing 
/2 contribution and the increasing pomeron contribution, 
whose trajectory has an intercept of a(0) ~ 1-l-e. Although 
the JML model reproduces the general W dependence of 
the longitudinal exclusive p° cross section relatively well, 
it drops as a function of Q^ faster than the data and agrees 
only up to Q2 « 4.10 GeV^. 

We now turn to the GPD approaches. The dashed line 
shows the result of the GK model, while the thin solid 
line shows the result of the VGG model. We see that they 
give a good description of the high and intermediate W 
region, down to VF « 5 GeV. This result was already ob- 
served by the HERMES collaboration PD]. At high W the 
slow rise of the cross section is due to the gluon and sea 
contributions, while the valence quarks contribute only 
at small W (this decomposition is shown in fig. [55] when 
we discuss the transverse cross section). We see a signif- 
icant disagreement between the GK and VGG models at 
intermediate W, which can be clearly explained by the 
fact that, as was mentioned in section \3l2\ the GK model 
takes the interference between the two handbag diagrams 
of fig. [M]into account, while the VGG model sums them 
incoherently. This interference is of course maximal at in- 
termediate VF's where the gluon handbag diagram starts 
to become significant, while the valence part of the quark 
handbag diagram is still significant. The data don't par- 
ticularly favor GK over VGG but it is clear that, on purely 
theoretical grounds, the GK model is more correct. It is 
remarkable that, except in this intermediate W region, i.e. 



in the high- and low-IF regions, the GK and VGG mod- 
els are in close agreement. The fact that two independent 
groups with different numerical methods and approxima- 
tions tend to agree gives some relative confidence in the 
calculations. 

At lower W values, where the new CLAS data lie, it 
is striking that both the GK and VGG models fail to re- 
produce the data. This discrepancy can reach an order 
of magnitude at the lowest W values. The trend of these 
particular GPD calculations is to decrease as W decreases, 
whereas the data increase. In the VGG and GK calcula- 
tions, these trends can be understood as follows: GPDs 
are approximately proportional to the forward quark den- 
sities q{x). This relation is not so direct since the quark 
densities are, in the double distributions ansatz of ref. [46^ , 
convoluted with a meson distribution amplitude but, still, 
the main trends remain. Then, as x increases (i.e. W de- 
creases), GPDs tend to go to since q{x) « (1 — x)^ for x 
close to 1. There might be, according to the scale, a slight 
local increase or bump around x ~ 0.3, due to the va- 
lence contribution, which is indeed clearly apparent in the 
VGG calculation shown in fig.[l5] However, this variation 
cannot explain an increase of an order of magnitude. 

The conclusion on the GPD approach is then two-fold: 

— The handbag is not at all the dominant mechanism 
in the low W valence region and higher twists or so 
far uncontrollable non-perturbative effects obscure the 
handbag mechanism. If so, one has to explain why the 
(power-corrected) handbag mechanism works in the 
high/intermediate W (i.e. low x) domain and, quite 
abruptly, fails in the valence region. Higher twist can 
certainly depend on energy, but such a strong variation 
with W is certainly puzzling. Also, the explanation 
might simply be of a kinematic nature. As shown in 
fig. 1141 the minimum value of | i | increases significantly 
with decreasing W. For instance, tmin ~ 1-6 GeV^ for 
the (0.64 < XB < 0.70, 5.10 < Q^ < 5.60 GcV^) bin 
while tmin « 0.1 GeV^ for the (0.16 < xb < 0.22, 
1.60 < Q2 < 1.90 GeV^) bin. In the handbag formal- 
ism, higher twists grow with t and this purely kine- 
matic effect provides a natural source for them. How- 
ever, more than absolute values, the ratio -^ should be 
relevant, and for the largest tmin values, one actually 
finds ^^5^ = jj^, i.e. of the order of 30%. More gener- 
ally, the largest tmin values correspond to the largest 
Q^ values but, since Q^ increases faster than tmin in 
our kinematics, this actually makes the ratio ^ more 

favorable as Q^ increases. 

— Or the handbag mechanism, which succesfuUy describes 
the region of intermediate and high W, is indeed at 
work in the valence region but the way the GPDs are 
modeled by the VGG and GK groups is incomplete, 
with a significant and fundamental contribution miss- 
ing, or incorrect. 

Let us note at this stage that, on general grounds, at 
large xb , the situation is somewhat more complex than at 
small Xb- Both real and imaginary parts contribute (the 
so-called ERBL and DGLAP regions), the skewness of the 
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Fig. 25. World data for the reduced cross sections 7£p -^ ppl as a function of W for constant Q^ bins, in units of fiharn. The 
lowest cross section point in the 2.80 GeV < Q < 3.10 GeV bin (from CORNELL) corresponds to the low Rp (=0.38) point 
in fig. [221 and might be unreliable. The dashed curve shows the result of the GK calculation and the thin solid curve shows 
the result of the VGG calculation. Both calculations are based on Double Distributions as proposed in ref. ^S] for the GPD 
parametrizations and incorporate higher twist effects through k± dependence. They differ essentially in summing coherently or 
not the gluon and the quark exchange handbag contributions (see fig. I24|l . The thick solid curve is the VGG calculation with 
the addition of the D-term inspired contribution. The dot-dashed curve shows the results of the Regge JML calculation. The 
4.2 GeV CLAS, CORNELL, HERMES and E665 data are respectively from refs. [21], [19], [20] and i33j. 



GPDs is substantial and non-perturbative effects are ex- 
pected to play a strong role in determining the behavior 
of the GPDs near x ^ ^. Therefore, it is not so clear 
whether the higher twist corrections, which are already 
substantial at low xb through the k± dependence, have 
the same character at large xg and if much can be con- 
cluded about large xb from the good GPD description of 
the small xb data. 

Therefore, with utmost caution, we quote the sugges- 
tion of ref. [in] to add a new (strong) component to the 
standard VGG GPD parametrization, in the form of a D- 
term inspired ansatz, to reconcile the handbag approach 
with the data. We recall that the D-term was originally in- 
troduced by Polyakov and Weiss [50j in order to complete 
the Double Distribution representation of GPDs, so as to 
satisfy the polynomiality rule, and that it could be inter- 
preted as the contribution to the GPDs of the exchange of 
a a meson in the t-channel. In ref. [15] , the i-dependence of 
the D-tevm was modified (making it, effectively, no longer 
a D-term, properly speaking) and renormalized. One of 
the motivations for this new term was to extend the con- 
cept of qq components, or t-channel meson exchange, in 
GPDs, in a spirit similar to the JML model that explains 



the strong rise of the cross section as W decreases by t- 
channel a and /2 meson exchange processes. The thick 
solid line of fig. [551 shows the result of the introduction 
of this new contribution, added coherently to the stan- 
dard VGG double distribution parametrization, with its 
normalization adjusted to the data. 

We insist that this extra contribution is a speculation, 
which however does have the merit of providing numerical 
estimates of the cross section. Several alternative expla- 
nations should also be in order. GPDs can obviously be 
parametrized differently than in VGG and in GK. It was 
shown for instance that the spectator model of Hwang and 
Miiller with an overlap representation for the modeling 
of the GPDs [51] produces naturally an enhancement at 
large x compared to the VGG model that should produce 
quark exchange cross sections dropping with increasing 
W. Also, NLO QCD corrections might be more sizeable 
in this region (see ref. jS^)- Let us finally mention that, 
in the framework of the VGG model, the Feynman mech- 
anism (or overlap diagram, see ref. [15j for instance) was 
calculated but, although it has the right W dependence, 
since it is a real contribution that lives in the | a; |< ^ 
region and therefore dominates at large xb , its numerical 
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contribution, which does not rely on any extra parameter, 
is barely significant. 

In summary, the only conclusion that we allow our- 
selves to reach at this stage is that the two popular GK 
and VGG models that provide numerical estimations of 
the ep — *■ epp" cross sections and which describe well these 
data at large W values, fail to describe the present large 
Xb CLAS data (the normalisation as well as the W depen- 
dence). This fact does not imply however that the hand- 
bag mechanism is not at work in this latter regime as this 
might simply be an artifact of the current double distri- 
bution parametrization by these two groups. Wc reiterate 
that the necessary consistent treatment of GPD modeling, 
QCD scale setting, and higher-twist effects is much more 
difficult at low W than at high W and renders conclusions 
more difficult. 

Fig. 1261 shows the longitudinal diEev&it\a\ cross section 
daL/dt as a function of t. The model of ref. [49], inspired 
partly by Regge theory, naturally explains the decrease of 
the i-slope as xb gets larger. So does the spectator model 
of ref. [51] , which quotes a decrease of the t slope from h « 



Let us mention here that because of the statistics of 
the data, which in general fall quite rapidly with t, we have 
not been able to extract reliably the ratio Rp = ^^ as a 

function oft and that, in a given {xB^ Q^) bin, the same Rp 
value has been applied to the data over the whole t range. 
In other words, these differential cross sections might need 
to be corrected if the ratio Rp depends significantly on t. 

The dash-dotted curves in fig.l^show the JML model. 
In this approach, the natural t dependence given by the 
Regge formula s"*^*^ is too sharp, as for the H and E GPD 
case. Thus, i-dependent form factors have to be introduced 
at the electromagnetic vertices of the diagrams of fig. [23l 
according to the procedure and motivations of ref. |42j . 

Fig. [27] shows the Q^ dependence of the longitudinal 



p cross section cr^ for different xb values. When more 
than three points are present in an xb bin, we fit the Q^ 
dependence to the function ^^, with the extracted val- 
ues n displayed in the figure. We recall that the handbag 
formalism predicts a value of 3 for n at asymptotically 
large Q^ values. However, a smaller coefficient, i.e. a flat- 
3.5 GeV'^ at xb=0.2 and Q^—2 GeV'^ to b ~ 1.5 GeV'^ at ter Q^ dependence, is expected at these low Q^ values 



XB=0.6 and Q^=5 GeV^. 



due to preasymptotic (fcj^) effects. The thick solid (blue) 



25 




0.22 <x„< 0.28 



0.28 <x„< 0.34 





,„2 -2.06±1.67 
OlO<Q ) 



0.40 < x„< 0.46 



0.46 <x„< 0.52 





2 -2.77 ±1.34 

aL"(Q ) 



0.58 <x„< 0.64 



0.64 <x„< 0.70 



-4i: 



i CLAS (5.754 GeV) o 



Q^ (GeV^) 



Q^ (GeV^) 



Q^ (GeV^) 



Fig. 27. Cross sections ctl for the reaction 7*p -^ pp^^ as a function of Q^ for different bins in xb, in units of /ibarn. Tlie solid 
blue curve is the result of the VGG calculation with the Generalized D-term including the k± correction. The dashed curve 
shows the leading twist (i.e without the k± correction) VGG handbag calculation (with the Generalized D-term) . The solid red 
curve is a fit to the data using a ^i- function. 



curves m fig. [27] show the results of the VGG calculation 
including the extra aforementioned term. The magnitude 
and shape of the data are reasonably reproduced. The k± 
effects in the calculation flatten the Q^ slope of the cross 
section. For comparison, the asymptotic result, i.e. with- 
out k± effects, is shown as the dashed (black) curve in 
fig-HZl One sees that its normalization is of course higher. 
Indeed, k± effects reduce the cross section by a factor 2 to 
5 depending on xb and Q^. Also, the Q^ dependence of 
the asymptotic result is steeper, precisely i. At asymp- 
totically large Q^ values, the two calculations, i.e. with 
and without k± effects, are expected to agree. 

To complete the interpretation of our data, we finally 
turn to the transverse part of the cross section. We show 
in fig. [28] the W dependence of the ^^p -^ ppj^ cross sec- 
tion for different Q^ bins. The JML model (dot-dashed 
curve) once again reproduces the general shape of the W 
dependence of this cross section. Though, quantitatively, 
it seems to overestimate by w 30% the GLAS data at low 
Q^, where it is expected to be the most valid. The agree- 
ment for the longitudinal part of the cross section was 
much better for this kinematical region. 

The transverse part of the cross section doesn't lend 
itself straighforwardly to a GPD interpretation since it 
is higher twist. However, the GK group, taking into ac- 
count fcj^ effects, has been able to extend its analysis of 



the longitudinal cross section to the transverse case |53j . 
They showed that retaining the quark transverse momenta 
regularizes the infrared singularities occuring in the trans- 
verse process. The dashed curve in fig. [28] shows the re- 
sult of this calculation. The high energy E665 data, and to 
some extent, the HERMES data are well reproduced, thus 
comfirming the approach. However, as for the longitudinal 
cross section, the low W GLAS data are completely un- 
derestimated. The VGG model has not yet been extended 
to the transverse case but it is clearly expected that the 
addition of the Generalized D-term to the transverse pro- 
cess produces the same effect as for the longitudinal one 
and might explain the rise of the cross section at low W . 
Finally, fig. [29] shows the transverse differential cross 
section dax/dt compared to the JML model. One finds 
that the JML model tends to overestimate the experimen- 
tal cross sections at the « 30% level, especially at large t 
values. 



4 Summary 



Using the GLAS detector at JLab, we have collected the 
largest ever set of data for the ep — > e'pp^ reaction in the 
valence region. We have presented the Q^ and xb (and 
W) dependences of the total, longitudinal and transverse 
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cross sections, as well as of the differential cross section in 
t. 

The unique features that we have observed are: 

— The W dependence of our data shows a clear decrease 
of the cross sections with increasing W, in contrast 
to the higher W data (HERMES, HI, ZEUS, E665), 
which show cross sections that tend to be flat or slowly 
rising with W, 

— The t dependence exhibits a varying slope with energy: 
the slope increases as xb decreases. In particular, our 
t dependences are almost flat at our largest xb values, 

— The cross sections decrease with Q^ as approximatively 
i, i.e. in a flatter way than what is predicted by the 
asymptotic handbag diagram. 

These data and features can be interpreted in two 
ways: 

— In terms of hadronic degrees of freedom, i.e. meson tra- 
jectory exchanges in the i-channel, following the JML 
Regge model. In order to reproduce the rather flat t 
dependences varying with Q^ and xb , electromagnetic 
form factors varying with Q^ and t are necessary. For 
the longitudinal part of the cross section, good agree- 
ment with the data is found up to Q^ « 4.10 GeV^. For 
the transverse part of the cross section, there seems to 
be an overestimation (by « 30%) of the cross section. 

— In terms of partonic degrees of freedom, i.e. quark 
handbag diagrams and Generalized Parton Distribu- 
tions. However, the GK and VGG calculations cannot 
provide the right W dependence of the cross section. 
This does not imply that GPDs cannot be accessed 
through exclusive p° electroproduction in the valence 
region but that possibly the way double distributions 
are modeled or the hard scattering amplitude is cal- 
culated in these two particular approaches should be 
modified or revisited. We stress that in exclusive me- 
son electroproduction the GPD modeling problem is 
convoluted with other issues such as the treatment of 
the QCD scale setting, higher-twist effects, the meson 
distribution amplitude, etc. , rendering conclusions dif- 
ficult. The present data will provide important input to 
improve our understanding of these fundamental QCD 
issues. 
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Fig. 28. World data for the reduced cross sections 'y^p — > pp^ as a function of W for constant Q^ bins, in units of fiharn. The 
dashed curve shows the resuh of the GK calculation. The solid curve shows the contribution of the valence part of the quark 
exchange handbag diagram (fig. [53] left) while the dotted curve shows the sum of the sea quarks part of the quark exchange 
handbag (fig. [24] left) and of the gluon exchange (fig. [24] right) contributions. The dot-dashed curve shows the results of the 
Regge JML calculation. The 4.2 GeV CLAS, CORNELL, HERMES and E665 data are respectively from refs. [21], [l9], [20] 
and [33l. 
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Fig. 29. Cross section dar/dt (in /ib/GeV^) for all bins in {Q^,xb) as a function off (in GeV^). The dash-dotted curve is the 
result of the JML model. 



